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tracked  vehicle  with  dual  inertial  mass  shakers  mounted  on  top  was  used  for  the  mobile 
source.  The  source  accurately  transmitted  the  shaker  signal  into  the  ground,  and  its 
mobility  made  it  a  practical  choice  for  field  operations.  It  excited  Rayleigh  waves,  as 
desired,  but  also  generated  undesirable  P-waves  and  was  not  found  to  be  directional.  It 
proved  incapable  of  finding  a  target. 

Improvements,  such  as  a  deploying  an  array  of  mobile  sources  and  a  stronger 
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I.  INTRODUCTION 


A.  MOTIVATION 

Buried  mines  hinder  operations  for  the  US  Armed  Forces.  Mines  do  so  by  three 
main  methods:  1)  inflicting  casualties  on  soldiers,  Marines,  and  sailors;  2)  damaging  or 
destroying  equipment;  and  3)  changing  battle  plans  due  to  the  perceived  threat  of  mines. 
They  have  proven  to  be  an  inexpensive,  effective  means  of  fighting  for  groups  on  the 
defensive.  The  threat  of  mines  is  especially  relevant  now  due  to  the  new  role  of  the  US 
Navy  and  Marine  Corps,  which  is  concentrated  on  going  out  and  finding  the  enemy. 
Offensive  forces  which  constantly  push  forward,  like  the  US  Armed  Forces,  provide  the 
perfect  target  for  land  mines.  177  (13  percent)  of  the  US  casualties  in  the  Gulf  War  of 
1991  were  due  to  Iraqi  landmines  and/or  UXO  (unexploded  ordnance)  (Stohl,  2002). 

Mines  also  are  a  large  problem  for  innocent  civilians  in  countries  which  have 
hosted  major  conflicts  over  the  last  half  century,  such  as  France,  Vietnam,  Cambodia, 
Yugoslavia,  Angola  and  Afghanistan.  Annually,  in  Angola  alone,  there  are 
approximately  6,000  new  amputees  due  to  landmines  (VVAF,  1999).  Worldwide, 
approximately  70  civilians  are  killed  every  day  due  to  landmines  that  were  left  over  from 
conflicts  (Sheetz,  2000). 

B.  CURRENT  METHODS  OF  MINE  DETECTION 

The  most  common  method  for  finding  land  mines  is  to  use  a  metal  detector.  This 
is  both  dangerous  and  slow.  Metal  detectors  require  the  operator  to  be  within  l-2m  of  the 
mine,  which  puts  the  operator  within  the  “kill  radius”  of  most  landmines.  Furthermore, 
metal  detectors  are  inaccurate.  They  give  a  high  number  of  false  “hits”  due  to  iron  that 
may  occur  naturally  in  the  ground.  Also,  they  are  useless  against  landmines  made  of 
something  besides  metal,  such  as  plastic  or  vulcanized  rubber. 

Dogs  can  also  be  used  to  find  mines,  but  they  must  be  even  closer  to  the  mine 

than  a  metal  detector.  The  obvious  problem,  from  a  strictly  military  viewpoint,  is 
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inefficiency.  Mine  hunting  dogs  require  time,  money,  and  effort  to  train.  It  simply 
would  not  be  worth  it  to  invest  that  much  in  training  a  dog  that  will  just  get  killed  very 
soon  after  becoming  operationally  useful.  They  get  tired  and  lose  effectiveness  after 
about  an  hour  or  two.  Also,  most  trained  sniffing  dogs  can  only  recognize  one  or  two 
types  of  explosives,  further  limiting  their  usefulness. 

Another  technology  under  testing  is  commercially  known  as  “Mongoose”.  It 
consists  of  giant  rope  net  grid  (about  half  of  a  football-field  size),  with  a  small,  shaped 
charge  at  each  junction  on  the  grid.  This  net  is  rocket-fired  over  a  suspected  minefield;  it 
spreads  out,  and  when  it  hits  the  ground,  the  shaped  charges  fire  into  the  ground, 
detonating  any  landmines  in  the  immediate  surrounding  area.  While  effective  if  shot  at 
specific  area  known  to  have  a  high  concentration  of  mines,  it  can  also  be  very  inefficient. 
Each  shot  is  expensive,  and  more  importantly,  it  does  not  find  mines,  it  just  detonates 
them  if  they  are  nearby.  This  could  lead  to  a  situation  where  dozens  of  nets  are  just  shot 
all  over  the  place,  in  hopes  of  finding  mines,  wasting  money  and  time.  If  used  in 
combination  with  a  mine  detection  technology,  the  Mongoose  system  could  be  effective. 

Ground  penetrating  RADAR  (GPR)  is  another  mine-detection  technology  under 
consideration.  It  is  based  on  sending  radio  waves  into  the  ground,  and  looking  for 
returns,  just  like  air  based  RADAR.  However,  GPR  is  very  limited.  Typical  operating 
frequencies  are  around  300  MHz,  which  results  in  very  short  detection  ranges  due  to  high 
attenuation  (Dolphin,  1997).  Also,  many  soils  are  electrically  conductive,  which  sharply 
decreases  the  range  of  GPR. 

C.  OBJECTIVES 


SriJirtH  l'  iiln  I'lH 


Figure  1 .  Seismic  SONAR  Concept 
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The  concept  of  seismic  SONAR  is  shown  above.  It  consists  of:  sending  energy 
into  the  ground,  which  propagates  as  seismics  waves;  and  listening  for  an  “echo”  signal 
from  the  target. 

The  goal  of  this  investigation  was  to  develop  and  test  a  mobile  source  for  a 
seismic  SONAR  mine  detection  system. 

In  meeting  the  goal,  it  was  necessary  to  fulfill  the  following  requirements: 

-  build  or  acquire  a  tracked  vehicle  that  a  vibration  shaker  system  could  be  mounted  on; 

-  test  the  rigidity  of  the  vehicle  to  see  if  it  was  capable  of  transferring  a  clear  signal  to  the 
ground; 

-  measure  the  received  signal  at  various  distances  that  could  be  transferred  into  the 
ground  through  the  vehicle; 

-  establish  a  beam  pattern; 

-  determine  the  types  of  waves  which  were  propagating; 

-  find  a  target. 


3 


THIS  PAGE  LEFT  INTENTIONALLY  BLANK 


4 


II.  SEISMIC  WAVE  PROPERTIES 


Seismic  waves  travel  through  the  ground.  They  place  a  strain  on  the  ground, 
which  acts  as  an  elastic  medium.  The  assumption  that  the  ground  acts  elastically  is  valid 
because  of  the  low  level  of  ground  deformation  which  occurs  (Geophysics,  St. 
Andrew’s).  In  an  earthquake,  the  first  wave  to  arrive  is  the  P-wave,  or  compressional 
wave.  If  there  is  a  sub-surface  reflection  of  the  P-wave,  then  that  reflection  is  next  to 
arrive.  After  that,  the  s-wave,  an  elastic  shear  wave  arrives.  Last  to  arrive  are  the  elastic 
surface  waves  (Seismo  Lab,  Cal  Tech). 

A.  BODY  WAVES 

There  are  two  types  of  body  waves:  compressional  (P-waves);  and  shear  (S- 

waves). 


1.  Compressional  Waves 


Compressional  waves  propagate  longitudinally,  i.e,  the  motion  of  the  ground  is  in 
the  same  direction  as  the  wave  propagation.  Particle  motion  is  accordian-like,  i.e., 
neighboring  particles  alternate  between  moving  closer  to  each  other  and  moving  further 
apart.  They  are  called  “P-waves”  because  they  arrive  first  (primary). 


(a)  P-wave 


■  compressions 


WMmmm 


undisturbed  medium 


—  dilations 


Figure  2.  Propagation  of  a  Compressional  Wave  (Brown  et.  al,  1981) 
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The  compressional  wave  propagates  with  speed 


k  +  4/j/3 


P 


(2.1) 


Vp  is  the  speed,  k  is  the  bulk  modulus  of  the  ground,  p  is  the  shear  modulus  of 
the  ground,  and  p  is  the  density  of  the  ground. 

P-waves  are  not  desired  in  seismic  mine  detection  because  they  are  not  confined 
near  the  surface  of  the  ground,  and  can  produce  multiple,  unwanted  reflections  from 
bottom  layer  boundaries  which  are  deeper  than  the  mine  is  buried.  This  results  in  many 
reflected  rays  arriving  at  sensors.  These  multiple,  reflected  ray  paths  make  it  difficult  to 
differentiate  between  the  initial  drive  signal,  the  bottom  bounce  reflection  signals,  and  the 
return  signal  from  the  target.  An  illustration  of  bottom  bounce  is  shown  below  in  Figure 
3. 


Top  Surface  Boundary 


Figure  3.  Bottom  bounce  reflection  of  a  P-Wave. 


2.  Shear  Waves 

Shear  waves  propagate  by  shear  strain,  with  the  ground  displacement 
perpendicular  to  the  direction  of  propagation  of  the  wave.  They  can  be  polarized  eith 
horizontally  (SH)  or  vertically  (SV).  They  are  called  “S-waves”  because  they  arrive 
second,  after  the  P-wave. 
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Figure  4.  Propagation  of  a  Shear  Wave  (Brown  et.  al,  1981) 


The  velocity  of  S-waves  is: 

V,=  p  (2.2) 

\P 

Vs  is  the  velocity,  p  is  the  shear  modulus  of  the  ground,  and  p  is  the  mass  density 
of  the  ground. 

S-waves  are  also  undesirable  for  mine  detection  because,  like  P-waves,  they  fill 
up  the  entire  homogeneous  half-space,  and  thus  can  produce  multiple  reflections. 

Both  types  of  body  waves  spread  proportional  to  1/r  within  the  ground,  where  r  is 
distance  from  the  source.  Along  the  surface,  they  spread  proportional  to  1/r2  (Gaghan, 
1998). 


B.  SURFACE  WAVES 


1.  Rayleigh  Waves 

Rayleigh  waves  travel  relatively  near  the  surface  of  the  ground.  The  particle 
motion  is  in  an  elliptical  orbit  in  a  vertical  plane  perpendicular  to  the  direction  of  wave 
propagation  (Geophysics  Dept.,  St.  Andrews).  The  vertical  and  horizontal  components  of 
particle  motion  are  90°  out  of  phase.  Near  the  upper  surface  boundary  (within  0.2 
wavelengths),  the  motion  is  retrograde,  i.e.,  the  particles  move  counter-clockwise  if  the 
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wave  is  propagating  from  left  to  right  (Fitzpatrick,  1998).  This  retrograde  elliptical 
motion  is  shown  by  the  darkened  ellipse  in  Figure  5.  At  depths  greater  than  0.2 
wavelengths,  the  particle  motion  switches  to  prograde,  with  the  particles  moving 
clockwise  as  the  wave  propagates  left  to  right  (Fitzpatrick,  1998),  as  shown  in  Figure  6. 


Rayleigh  Wave 


Figure  5.  Propagation  of  a  Rayleigh  Wave  (Brown,  et.  al,  1981) 


PROGRADE 


RETROGRADE 


Vertical  Signal  Leading  Radial 
Signal  by  90  degrees 


Radial  Signal  Leading  Vertical 
Signal  by  90  degrees 


Figure  6.  Prograde  and  Retrograde  Rayleigh  waves. 


Rayleigh  waves  propagate  with  velocity: 

cr  =  csk 


(2.3) 
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where  cr  is  the  Rayleigh  wave  speed,  cs  is  the  shear  wave  speed,  and  k  is  a  dimensionless 
parameter,  given  by: 

K6  -8/r4  +(24-16r)w2  +(16 r2  -16)  =  0 

2 c s  /u  _  l- 2v  (2.4  and  2.5) 

cl  A  +  2/j  2-2v 

where  y  is  a  dimensionless  parameter,  ci  is  the  longitudinal  wavespeed,  X  and  p  are  the 
Lame  constants,  and  v  is  Poisson’s  Ratio  (Fitzpatrick,  1998).  For  most  solid/fluid 
boundaries  (including  sand/air)  v  is  approximately  0.25,  and  k  is  approximately  0.9 
(Fitzpatrick,  1998). 

Rayleigh  waves  spread  out  with  range  proportional  to  1/  Vr  .  Recalling  that  body 
waves  spread  with  1/r,  over  large  horizontal  distance,  this  difference  in  spreading  rates 
results  in  Rayleigh  waves  being  the  dominant  form  of  seismic  energy.  They  also  decay 
very  rapidly  with  depth,  carrying  almost  no  energy  deeper  than  2  wavelengths  (Seismo 
Lab,  Cal  Tech).  This  depth  decay  is  ideal  for  seismic  SONAR  because  it  means  that 
Rayleigh  waves  will  not  give  strong  reflections  off  of  bedrock,  etc.,  which  would 
complicate  the  interpretation  of  return  signals. 


2.  Love  Waves 

Love  waves  are  ducted,  horizontally-polarized  shear  waves.  Their  oscillations  are 
in  the  plane  of  the  surface  of  the  ground,  and  perpendicular  to  the  direction  of 
propagation  (see  Figure  7).  There  was  no  interest  with  Love  Waves  for  this  project 
because  they  vibrate  horizontally  with  respect  to  the  ground.  The  source  for  this  project 
vibrated  vertically  with  respect  to  the  ground,  hence  Love  waves  would  only  be 
minimally  excited,  if  at  all.  However,  they  still  could  be  used  for  mine  detection  with 
sources  that  vibrate  horizontally  with  respect  to  the  ground.  Their  main  drawback  is  that 
they  are  a  waveguide  mode  and  possess  multiple  cutoff  frequencies.  They  are  also 
dispersive. 
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Love  Wave 


Figure  7.  Love  Wave  Propagation,  left  to  right  (Brown,  et.  al,  1981) 


C.  SEISMIC  WAVES  FROM  A  VERTICAL  SOURCE 

Seismic  wave  interactions  at  a  fluid-solid  boundary  depend  on  the  boundary 
conditions  (moduli,  layering,  etc.).  These  conditions  determine  the  propagation  speed, 
dispersion,  depth  of  vibration  in  the  sediment,  and  height  of  vibration  in  the  air.  For  this 
research,  the  sand/air  boundary  resulted  in  negligible  energy  being  transferred  into  the 
air.  After  a  short  settling  distance,  the  resulting  wave  is  a  “smooth  band”  of  seismic 
energy  which  propagates  radially  and  vertically.  This  has  been  discussed  in  much  greater 
detail  in  previous  theses  (Gaghan,  Fitzpatrick,  etc.) 

Using  geophones,  it  is  possible  to  detect  the  ground  disturbance  caused  by  seismic 
waves.  Seismic  waves  also  scatter  from  strong  enough  inhomogeneties,  such  as  a  mine¬ 
like  target,  and  can  be  detected.  The  ability  to  detect  these  waves  is  the  basis  for  a 
seismic  SONAR  concept.  As  the  figure  below  shows,  ideally,  the  majority  of  energy  that 
the  geophones  detect  is  in  the  form  of  Rayleigh  waves.  A  major  assumption  of  Fig.  8 
below  is  that  of  a  semi-infinite  half-space,  with  no  bottom  interactions.  If  a  bottom  exists 
at  a  shallow  enough  depth  where  P-waves  are  still  relatively  strong,  bottom  bounce  of  P- 
waves  could  occur.  Bottom  bounce  is  when  a  wave  bounces  off  the  bottom,  goes  to  the 
top  surface,  reflects  off  that,  and  continues  this  pattern  of  bottom/top  reflection,  moving 
horizontally  all  the  while.  A  diagram  of  bottom  bounce  is  shown  above  in  Figure  3. 
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Circular  Footing 


- - 

Wave  Type 

Per  Cent  of 

Total  Energy 

Rayleigh 

67 

Shear 

26 

Compression 

7 

Figure  8.  Composite  of  all  waveforms  from  a  circular  footing,  in  a  homogeneous,  isotropic,  elastic  half¬ 
space  (Richart,  et.  al) 
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III.  PRIOR  RESEARCH 


A.  UNIVERSITY  OF  TEXAS 

Research  on  the  seismic  SONAR  was  initiated  by  Dr.  Thomas  Muir,  then  at  the 
Applied  Research  Laboratories  of  the  University  of  Texas  at  Austin  (ARL:UT)  in  the 
1990-1997  timeframe.  This  work  was  continued  at  NPS  by  Prof  Muir  when  he  assumed 
the  chair  professorship  of  Mine  Warfare  in  the  late  1990’ s.  The  research  group’s  goal 
was  to  see  if  mines  could  be  detected  on  a  beach,  for  eventual  applications  to  amphibious 
landings.  The  source  used  was  a  heavy  rectangular  plate  with  nails  sticking  out,  and 
driven  by  an  electro-mechanical  transducer.  They  then  set  up  seismometers  and 
attempted  to  find  a  large  metal  cylinder. 

While  out  in  the  field,  they  acquired  echo  ranging  data  for  subsequent  analysis. 
However,  when  they  returned  to  their  lab,  they  were  able  to  detect  the  target  after 
extensive  signal  processing.  The  signal  processing  included  coherent  averaging, 
background  subtraction,  and  vector  polarization  filtering  (Fitzpatrick,  1998). 

B.  NAVAL  POSTGRADUATE  SCHOOL 

1.  LT  Stewart,  USN 

The  first  research  at  NPS  was  done  by  LT  William  Stewart  in  1995  (Stewart, 
1995).  He  created  a  conical-shaped  plunger-type  source  using  a  loudspeaker  mounted 
just  above  the  ground,  and  tested  its  transmitted  signal  over  a  large  range  of  frequencies. 
His  testing  was  done  in  the  large  sand  tank  (an  above  ground  swimming  pool  filled  with 
sand)  in  the  basement  of  Spanagel  Hall.  Using  seismometers,  he  was  able  to  show  that 
his  source  could  put  a  relatively  strong  seismic  signal  into  the  ground.  However,  the  tank 
was  too  small  for  echo  ranging  experiments. 
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2.  LT  Gaghan,  USN 


Student  Frederick  Gaghan,  with  advisors  Steven  Baker  and  Thomas  Muir, 
continued  the  research  at  NPS  in  1998  (Gaghan,  1998).  His  focus  was  on  source 
development  and  discrete-mode  excitation. 

His  source  consisted  of  two  inertial  mass  shakers.  He  mounted  them  on  an 
aluminum  framework,  so  that  each  was  pointing  down  at  45°.  He  had  them  at  45° 
because  he  needed  two  degrees  of  excitational  freedom  to  excite  elliptical  Rayleigh 
waves.  He  was  successful  in  preferentially  exciting  Rayleigh  waves.  He  reported  that 
discrete-mode  excitation  of  Rayleigh  waves  had  high  potential,  but  a  better  source  would 
improve  the  preferential  Rayleigh  excitation. 

3.  LT  Fitzpatrick,  USN 

Student  Sean  Fitzpatrick,  with  advisors  Steven  Baker,  Thomas  Muir,  and 
Anthony  Healey  picked  up  where  LT  Gaghan  left  off  (Fitzpatrick,  1998).  His  focus  was 
on  improving  source  design  from  LT  Gaghan’ s  work.  He  ended  up  using  a  linear 
magnetic  force  actuator,  which  was  also  an  intertial  mass  shaker,  with  a  metal  plate 
attached  to  a  rod  that  pounded  on  the  ground.  He  made  the  source  water-tight  so  that  it 
might  have  applications  for  mines  buried  in  shallow  water,  such  as  the  surf  zone. 

Using  a  two  element  seismometer  array,  he  was  able  to  locate  71  -  291  kg  targets 
at  ranges  up  to  5m. 

4.  MAJ  Hall,  USMC 


Student  Patrick  Hall  focused  his  research  on  target  strength  (Hall,  1998).  He 
worked  together  with  LT  Fitzpatrick.  He  measured  the  reflectivity  of  targets  as  a 
function  of  their  mass  load.  He  found  that  target  strength  increased  with  increasing  target 
mass. 
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5.  CPT  Sheetz,  USA 


Student  Kraig  Sheetz,  with  advisors  Thomas  Muir  and  Steven  Baker,  continued 
the  seismic  SONAR  work  at  NPS  in  2000  (Sheetz,  2000).  He  developed  a  new  source 
and  was  able  to  locate  an  M-19,  20  lb,  anti-tank  mine.  The  source  used  was  fourteen 
inertial  mass  shakers,  mounted  into  seven  pairs,  laid  out  in  a  linear  array  and  buried  in  the 
sand.  The  7-element  array  provided  high  directionality,  and  he  was  able  to  focus  a  large 
amount  on  energy  on  the  target.  He  was  able  to  demonstrate  signal  to  noise  ratios  of  9 
dB  in  the  echo  from  the  anti-tank  mine.  While  not  the  most  practical  of  source  designs, 
the  size  of  a  seven  element  array  was  good  for  detection,  the  concept  of  mine  detection 
with  seismic  waves  received  a  large  boost  of  confidence. 

6.  LT  McClelland,  USN 

Student  Scott  McClelland  did  work  on  rolling  seismic  SONAR  source 
(McClelland,  2002).  His  source  consisted  of  a  two  shaker  mounted  in  a  manually-pushed 
rolling  cylinder.  He  was  able  to  detect  a  1000-lb  bomb  at  5  meters.  The  roller  could  only 
take  readings  when  the  shakers  were  directly  aligned  with  the  ground  (i.e.,  once  per 
cylinder  revolution),  and  thus  was  not  ideal. 
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IV.  EQUIPMENT 


A.  SOURCE  DESIGN 

It  is  essential  for  a  deployable  seismic  SONAR  system  to  have  a  mobile  source. 
Simply  stated,  it  is  impractical  to  bring  a  system  out  into  the  field  if  it  cannot  be  easily 
moved.  With  speed,  adaptability,  and  efficiency  being  vital  in  today’s  battlefield,  mobile 
source  development  is  key  in  advancing  the  field  of  seismic  SONAR. 

Many  ideas  were  considered,  including:  placing  the  source  in  a  manually-pushed 
roller;  placing  the  source  on  a  sled,  which  could  be  pushed  by  an  ATV;  and  placing  the 
source  on  top  of  a  tracked  vehicle.  After  careful  consideration  of  each,  it  was  decided 
that  mounting  a  source  on  top  of  a  tracked  vehicle  would  be  the  best.  A  tracked  vehicle 
offered  benefits  such  as:  durability;  and  the  ability  to  handle  almost  any  terrain;  compact 
size;  and  it  could  be  remotely  controlled,  keeping  operators  away  from  mines. 

Fortunately,  a  tracked  vehicle  was  already  available  for  use.  This  made  the  work 
considerably  easier  because  construction  of  such  a  vehicle  would  have  set  us  back  several 
months.  Professor  Richard  Harkins  had  obtained  a  remote  controlled  tracked  vehicle  for 
his  robotics  class,  which  he  allowed  us  to  use. 

The  robot  was  capable  moving  forward  and  backwards,  turning  left  and  right,  and 
could  be  controlled  remotely.  However,  all  of  the  testing  was  performed  with  the  robot 
stationary.  The  idea  was  that  the  robot  could  go  to  a  location,  send  out  seismic  waves  and 
look  for  mines,  then  move  to  another  location  and  look  for  mines  there.  The  robot  is 
capable  of  sending  out  enough  signals  in  a  short  period  of  time,  to  keep  the  speed  and 
efficiency  of  mine  detection  high. 
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Figure  8.  Robot  with  2  shakers  mounted.  The  LCD  screen  is  a  digital  compass.  Upper  right  corner  shows 
a  3-axis  geophone. 
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The  tracked  vehicle  used  was  27  cm  wide  (cog  to  cog),  40  cm  long,  and  14  cm 
high  (ground  to  top  of  chassis).  Each  cog  was  16  cm  in  diameter.  The  robot  mass  was  12 
kg  when  no  equipment  was  mounted  on  it. 

The  robot  moved  by  means  of  4  fiberglass  reinforced  plastic,  cogged  wheels,  with 
plastic  tracks  on  each  side.  The  body  of  the  robot  consisted  of  wiring,  electronics,  etc., 
with  aluminum  support  rods,  and  outer  aluminum  plating  to  act  as  a  shell. 

A  1  cm  thick,  14  cm  x  14  cm  square,  aluminum  plate  was  fastened  to  the  top  of 
the  chassis  of  the  robot,  over  the  center  of  the  front  axle.  Then,  two  Aura  Bass  Shakers 
(inertial  mass  shakers)  were  mounted  on  the  plate,  one  shaker  above  the  other.  The 
bottom  shaker  was  in  direct  contact  with  the  plate,  and  the  upper  speaker  was  mounted  on 
aluminum  rods  directly  above  the  lower  shaker.  The  shakers  were  wired  in  series,  and  so 
they  “thumped”  at  the  same  time  and  in  phase  with  each  other.  Each  shaker  was  1.2  kg, 
giving  a  total  robot  mass  of  14.4  kg. 

The  shakers  were  centered  on  the  axle  because  the  axle  was  the  most  rigid  part  of 
the  robot,  and  it  also  offered  the  most  direct  path  of  coupling  between  the  source  and  the 
cogs.  The  theoretical  assumption  was  that  the  robot  source  acts  as  a  two-element  array, 
with  each  front  cog  being  one  element  of  the  array.  In  order  to  get  a  two-element  array, 
the  cogs  and  the  source  had  to  be  coupled  as  well  as  possible.  The  rear  axle  of  the  robot 
was  only  confined  to  a  slot,  and  held  in  place  with  springs.  Therefore,  it  had  some 
“play”.  The  front  axle  of  the  robot  was  fixed  in  place,  and  so  the  “play”  was  negligible 
compared  to  the  rear.  Also,  by  confining  the  signal  to  the  front  of  the  robot,  it  kept  the 
footprint  length  down.  At  the  chosen  operating  frequency,  100  Hz,  the  Rayleigh 
wavelength  was  about  lm.  If  the  source  was  mounted  in  the  center  of  the  robot,  the 
footprint  would  have  been  the  length  of  the  robot,  or  about  half  a  wavelength.  By 
limiting  the  signal  transmission  to  the  width  of  the  robot  (the  front),  the  source  was 
comparatively  smaller  in  with  respect  to  the  produced  wavelength  (about  one-quarter 
wavelength). 

B.  TRAILER 

A  trailer,  designed  and  built  by  previous  students  working  on  seismic  SONAR, 
housed  all  of  the  equipment.  It  was  essentially  a  sheet- wood  box  on  top  of  a  metal  trailer. 
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The  headphone  pre-amplifier  was  an  RCA  SA-155  Integrated  Stereo  Amplifier. 
The  headphones  were  full  ear  cup  Sony  Noise-Canceling,  model  MDR-NC20.  The 
headphones  were  used  for  listening  to  signals  received  by  the  geophones. 

The  geophones  were  tri-axial  (x-,  y-,  and  z-direction).  This  was  achieved  by 
having  three  single-axis  geophones  in  a  single  housing.  A  representative  drawing  is 
shown  below.  Each  geophone  was  potted  in  rubber  and  was  waterproof.  Jay  Adeff 
constructed  them. 


Figure  10.  Equipment  Trailer. 
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Figure  11.  Drawing  of  a  single  geophone. 

A  calibration  sheet  is  shown  in  Appendix  B . 


The  Geophone  Panel  was  designed  by  LT  Stephen  Rumph,  and  constructed  by 
Samuel  Barone.  It  was  capable  of  receiving  input  from  up  to  18  geophones,  with  each 
geophone  transmitting  an  x-,  y-,  and  z-direction  signal.  The  geophones  could  be  grouped 
into  1  or  2  arrays.  The  panel  was  wired  so  that  each  geophone  was  in  parallel  with  the 
other  geophones.  Output  could  be  sent  anywhere  via  BNC  connection. 

The  power  generator  was  gas  powered  and  put  out  2000  W.  It  was  Honda,  model 

2000i. 

The  function  generator  created  a  drive  signal  for  the  transient  burst  measurements. 
It  was  a  Hewlett  Packard  model  3314A. 

The  frequency  filters  were  not  used,  as  one  did  not  even  work  and  the  other 
distorted  the  drive  signal  too  much. 

The  drive  signal  amplifiers  were  used  to  amplify  the  signal  out  of  the  function 
generator.  They  were  Carvin  model  DCM  1000,  with  1000  W  maximum  output. 

The  oscilloscopes  were  used  to  view  and  store  output  from  geophones  and 
accelerometers.  The  majority  of  the  transient  burst  measurements  were  done  with  a 
Tektronix  TDS  3014  Four  Channel  Color  Digital  Phosphor  Oscilloscope  (pictured).  For 
frequency  scan  response  measurements,  a  Stanford  Research  Systems  SR785  Signal 
Analyzer  (not  pictured)  was  used.  For  the  target  detection  measurements,  a  Hewlett 


21 


Packard  Infinium  Two  Channel  Oscilloscope  (not  pictured)  was  used.  It  was  used  instead 
of  the  Tektronix  because  it  had  less  internal  electrical  “noise”. 

When  on  the  beach,  the  trailer  was  towed  with  a  John  Deere  “Trail  Gator”. 
Typical  cars  were  not  capable  of  driving  through  the  sand,  so  the  “Gator”  was  necessary, 
and  it  performed  well. 
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V.  TESTING 


Testing  was  conducted  from  February  2003  -  May  2003.  For  all  of  February,  and 
March,  the  testing  was  done  in  the  sand-filled  outdoor  swimming  pool  located  in 
Spanagel  Hall  Room  025.  In  April  and  May,  testing  was  done:  in  the  sand-filled  tank;  on 
the  NPS  softball  field;  and  at  Del  Monte  Beach. 


A.  TANK  ROOM 

1.  Source  Choice 

The  first  task  was  to  select  a  shaker.  There  were  two  options:  a  “Buttkicker  2”, 
made  by  The  Guitammer  Company;  or  dual  “Bass  Shakers”,  wired  in  series,  made  by 
Aura  (pictured  in  “Source  Design”  chapter). 

For  the  purposes  of  testing  the  shakers,  and  for  many  subsequent  tests,  a  transient 
burst  signal  was  used.  This  signal  was  a  single  sine  wave,  at  a  given  frequency  and 
voltage,  generated  once  per  second.  A  qualitative  representation  of  the  drive  signal  used 
for  testing  is  shown  in  Figure  12. 


Him.  H 

Figure  12.  Representation  of  drive  signal  used  in  transient  burst  measurements. 
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As  can  be  seen,  there  is  a  single  sine  wave  burst  at  the  beginning  of  each  second. 

The  first  source  test  was  on  the  Buttkicker  2.  The  Buttkicker  2  is  intended  for  the 
entertainment  industry,  to  shake  the  floor  in  order  to  create  the  sensation  of  powerful  bass 
in  music.  It  is  an  electro-magnetic  transducer  that  moves  a  3.5  lb.  inertial  mass  (a 
magnet)  up  and  down.  The  mass  is  suspended  magnetically,  with  no  mechanical 
connections.  It  has  a  wide  frequency  range,  from  approximately  5  Hz  to  200  Hz 
(Chiarealla,  2001). 

Figure  13  shows  the  acceleration  response  of  the  robot  for  when  a  Buttkicker  2 
was  mounted  on  the  on  top  of  the  chassis  of  the  robot.  Two  PCB  Piezotronics  shear 
accelerometers,  model  J353B03,  were  used  to  record  motion  caused  by  the  Buttkicker  2. 
Specifications  for  the  accelerometers  are  shown  in  Appendix  C.  One  accelerometer  was 
mounted  on  the  chassis  directly  above  the  front  axle  of  the  robot,  and  the  other 
accelerometer  was  buried  directly  beneath  the  left  front  cog,  approximately  3  inches 
deep.  The  Buttkicker  was  driven  at  63.6  V  (peak  to  peak)  with  a  sine  wave  burst  (one  per 
second)  at  100  Hz.  The  drive  voltage  was  chosen  so  that  acceleration  would  be 
approximately  lg.  It  was  desired  to  keep  robot  chassis  acceleration  due  to  the  shakers 
under  1  g  so  that  the  robot  would  not  jump  off  the  ground.  Any  jumping  off  the  ground 
would  cause  non-linearities  by  sending  an  un-controlled  wave  into  the  ground,  and  would 
make  signal  processing  more  difficult.  The  same  test  was  conducted  at  50  Hz  and  70  Hz, 
with  similar  results. 

The  y-axis  shows  acceleration  response,  in  terms  of  [g],  with  lg  =  9.8 1  m/s",  and 
the  x-axis  shows  time  in  [s] .  This  plot  is  for  a  single  sine  wave  drive  burst.  The 
averaging  function  of  the  Tektronix  oscilloscope  was  used.  Thirty-two  averaging  was 
chosen,  which  meant  that  the  oscilloscope  displayed  an  average  of  the  last  32  signals  it 
received,  as  opposed  to  showing  each  signal  as  it  happened.  The  averaging  function 
proved  useful  in  clearing  up  random  noise  and  interference.  Thirty-two  averaging  was 
used  on  all  testing  involving  using  the  Tektronix  oscilloscope.  The  Tektronix  was  used 
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for  measuring  all  transient  burst  measurements,  unless  otherwise  noted.  The  response  of 
the  plate  (black)  is  the  signal  with  larger  amplitude.  The  sand  response  is  shown  in  gray. 
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Figure  13.  Robot  and  sand  response  with  a  ButtKicker  2  as  a  source. 


Figure  13  shows  that  the  Buttkicker  2  failed  to  deliver  a  “clean”  signal.  There 
was  a  lot  of  signal  fluctuation  on  the  robot  around  0.02  seconds.  These  large  signal 
fluctuations  indicate  non-linearities,  which  would  complicate  signal  processing.  Also, 
the  signal  from  the  Buttkicker  2  coupled  poorly  to  the  sand.  There  was  approximately  a 
95%  drop-off  in  acceleration  amplitude.  This  plot  indicates  that  the  Buttkicker  2  is  not  a 
good  source  for  testing,  since  the  goal  was  a  relatively  clean,  regular  signal. 

The  other  choice  was  the  Aura  Bass  Shaker.  Manufacturer’s  specifications  are  in 
Appendix  A.  It  is  an  electro-magnetic  transducer.  Similar  to  the  Buttkicker  2,  it  moves 
a  magnet  up  and  down  to  generate  vibration. 

Figure  14  below  shows  information  for  the  same  experimental  set-up  as  Figure  13 
above,  but  dual  Aura  Bass  Shakers  were  mounted  on  top  of  the  chassis  instead  of  a 
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Buttkicker  2.  The  drive  signal  was  lower  than  for  the  Buttkicker  2  test  (41.5  V  peak  to 
peak  vs.  63.6  V),  but  acceleration  was  kept  the  same  at  approximately  lg.  Robot 
acceleration  response  is  shown  in  black,  and  sand  acceleration  response  is  in  gray. 


Figure  14.  Robot  and  Sand  acceleration  response  with  dual  Aura  Bass  Shakers  as  source. 


The  x-axis  shows  time  in  [s],  and  the  y-axis  shows  acceleration  response  in  [g] . 
As  can  been  seen  from  the  plot,  the  dual  Aura  Bass  Shaker  set-up  delivered  a  much 
cleaner  signal  than  the  Buttkicker  2.  There  was  no  significant  distortion  anytime  during 
the  drive  signal.  Also,  the  signal  transmission  to  the  sand  was  stronger  with  Aura  Bass 
Shakers.  This  plot  shows  approximately  a  70%  drop-off  in  acceleration  amplitude  from 
the  robot  to  the  sand,  while  the  Buttkicker  2  showed  approximately  95%  drop-off.  The 
gradual  decrease  in  force  amplitude  over  time  in  this  plot  indicates  that  there  was  some 
resonant  ringing,  which  could  not  be  avoided. 

Before  officially  determining  that  Aura  Bass  Shakers  would  be  used  as  the  source, 
testing  was  done  for  further  confirmation. 
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The  first  test  involved  seeing  how  much  signal  attenuation  occurred  in  the  robot 
itself  with  Aura  Bass  Shakers  as  source.  This  was  accomplished  by  running  two  tests, 
both  of  which  had  one  accelerometer  mounted  on  top  of  the  top  shaker,  and  one 
accelerometer  buried  approximately  3  inches  beneath  the  surface  of  the  sand.  The  first 
test  had  the  shakers  mounted  on  the  robot,  and  the  second  test  had  the  shakers  directly  on 
the  ground.  On  the  “with-robot”  experiment,  the  sand  accelerometer  was  directly  beneath 
the  left  front  cog.  On  the  “no-robot”  experiment,  the  sand  accelerometer  was  directly 
beneath  the  center  of  the  shakers.  The  results  are  shown  below. 
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Figure  15.  Sand  and  Shaker  acceleration  response  when  shaker  mounted  on  robot. 


Figure  15  shows  the  results  for  the  shakers  mounted  on  the  robot.  Time  in  [s]  is 
shown  on  the  x-axis,  and  acceleration  response  in  [g]  on  the  y-axis.  Speaker  acceleration 
response  is  in  black,  and  sand  acceleration  response  is  in  gray.  Drive  voltage  was  33.0  V, 
peak  to  peak.  Accurate  transmission  of  phase  and  shape  is  shown,  with  a  signal 
amplitude  drop-off  of  approximately  85%.  Resonant  ringing  was  present. 


27 


Figure  16  shows  force  response  curves  for  when  the  shakers  were  directly  on  the 
ground.  These  data  are  not  shown  on  the  same  plot  as  the  robot  curve  because  it  would 
be  too  much  to  show  on  one  plot.  The  data  would  be  unreadable. 
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Figure  16.  Shaker  and  Sand  acceleration  response  for  shaker  mounted  directly  on  ground. 


This  figure  shows  time  in  [s]  on  the  x-axis  and  acceleration  response  in  [g]  on  the 
y-axis.  Acceleration  response  is  shown  in  black,  and  sand  acceleration  response  is  shown 
in  gray.  Signal  transmission  was  very  good.  Shape  and  phase  were  maintained,  and  the 
drop-off  in  acceleration  amplitude  was  approximately  50%.  There  was  some  resonant 
ringing  from  the  shaker  case,  but  not  as  much  as  in  the  “with-robot”  experiment. 

Figure  17  shows  a  comparison  of  acceleration  response  of  the  robot  vs. 
acceleration  response  of  the  shaker.  Robot  acceleration  response  is  shown  in  gray,  and 
shaker  acceleration  response  is  shown  in  gray.  The  robot  accelerometer  was  mounted  on 
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the  bottom  of  the  front  axle.  Since  it  was  mounted  upside-down,  a  polarity  correction 
had  to  be  made,  which  was  just  multiplying  the  received  signal  by  (-1). 


All  Acccltrc  utters  100  Hi.  aa.C  V  Dr  t  F*nk  to  Peak.  m  flHxr.  averaging.  orlh  polarity 


Figure  17.  Shaker  and  robot  force  response  curves. 


Time  in  [s]  is  on  the  x-axis,  and  acceleration  response  in  [g]  is  on  the  y-axis. 
Almost  perfect  signal  transmission  was  shown  through  the  body  of  the  robot.  Phase  and 
shape  were  well  maintained,  with  approximately  a  15%  drop-off  in  force  amplitude.  The 
same  ringing  was  present  as  was  present  in  previous  experiments  involving  robot- 
mounted  shakers. 

These  last  three  plots  showed  that  the  robot  was  responsible  for  both  signal 
attenuation  and  some  resonant  ringing  not  present  when  the  shakers  were  directly  on  the 
ground.  The  third  plot  showed  almost  perfect  signal  transmission  from  the  shaker  into 
the  body  of  the  robot,  but  the  first  plot  showed  significant  decrease  of  the  signal  through 
the  robot  by  the  time  it  reached  the  sand.  That  indicated  that  the  robot-caused  losses 
occurred  in  between  the  body  and  the  sand.  The  only  objects  in  between  the  body  of  the 
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robot  and  the  sand  were  the  suspension  springs,  bearings,  and  the  cogs.  While  the  losses 
are  unfortunate,  they  are  unavoidable  since  suspension  is  needed  for  the  robot  to  be  an 
all-terrain  vehicle. 


2.  Source  Improvement 


The  Aura  Bass  Shakers  showed  a  relatively  clean,  repeatable  signal  transmission 
into  the  sand,  however  there  was  significant  power  loss.  One  of  the  major  requirements 
for  landmine  detection  was  a  strong  signal,  so  it  was  necessary  to  figure  out  a  way  to 
increase  the  signal  strength.  In  response  to  this  need,  the  dual  Aura  Bass  Shaker  set-up 
was  kept,  but  added  weight  on  top  of  to  the  chassis  of  the  robot.  After  trying  many 
different  weights  and  arrangements,  a  lead  brick  laid  on  top  of  the  chassis  of  the  robot 
proved  to  work  best.  The  brick  was  18.5  kg,  increasing  total  mass  of  the  robot  (robot, 
brick,  and  shakers)  to  33.3  kg. 

The  robot  accelerometer  was  mounted  on  the  bottom  of  the  front  axle,  and  the 
sand  accelerometer  was  buried  3  inches  beneath  the  left  front  cog.  Robot  response  is 
shown  in  black,  and  sand  response  is  shown  in  gray. 
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Figure  18.  Robot  and  Sand  force  response  for  dual  Aura  Bass  Shaker  mount  with  weight  added. 

As  Figure  18  clearly  shows,  there  was  a  much  improved  signal  transmission  into 
the  ground.  Signal  loss  was  reduced  to  approximately  50%,  down  from  the  70%  seen  in 
Figure  15.  A  higher  drive  voltage  was  used  than  for  an  “un- weighted”  robot,  but  robot 
acceleration  was  the  same.  The  end  result  was  that  more  force  was  put  into  the  ground 
when  a  mass  was  laid  on  the  chassis.  In  the  “un- weighted”  experiment,  using  F=ma,  the 
shakers  generated  1.05^  x9.81  mis2  / g  *13.4 kg  =  138A  of  force.  In  the  "with- weight" 
experiment,  the  shakers  generated  .80g  x  9.81m/ s2  /  g  x  33.3kg  =  26 UV  of  force.  If  the 
acceleration  was  taken  up  to  lg,  as  in  the  un-weighted  experiment,  326  N  could  have 
been  generated.  Clearly,  placing  weight  on  the  chassis  increased  force  generation  and 
improved  the  quality  of  the  source. 

Despite  the  amplitude  loss  and  resonant  ringing,  the  quality  of  signal  transferred 
into  the  ground  was  still  clean  and  repeatable,  and  thus  further  experimentation  with  the 
robot-mounted  source  was  worthwhile.  Dual  Aura  Bass  shakers  were  acceptable. 
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B.  SOURCE  PROPERTIES 


1.  Beam  Pattern 

The  first  task  in  determining  source  properties  was  directionality.  It  was  assumed 
the  robot  would  act  as  a  two-element  array,  with  each  front  cog  acting  as  one  element.  In 
order  to  test  this  assumption,  testing  was  done  at  the  NPS  softball  field.  The  trailer  was 
set  up  in  foul  territory  near  third  base,  and  the  robot  was  on  the  infield  near  second  base, 
approximately  25  meters  away  from  the  trailer.  The  infield  was  made  of  packed  dirt. 
One  geophone  was  placed  2m  away  from  the  robot,  and  another  geophone  was  4m  away 
(shown  below).  The  geophones  were  placed  on  the  surface  of  the  ground. 


Figure  19.  Photo  of  testing  set-up  on  softball  field. 


(not  to  scale) 
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Figure  20.  Schematic  of  testing  on  softball  field. 
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Ground  truthing  (sound  speed  measurements)  were  the  first  test  performed.  A 
single  impulse  was  generated  by  the  shakers,  and  each  geophone  received  the  signal.  The 
time  delay  between  when  each  geophone  received  a  signal  from  the  robot  was  recorded. 
The  distance  between  the  geophones  was  set  at  2m.  It  was  then  possible  to  calculate 
speed  of  sound  with  the  time  it  took  for  the  signal  to  get  from  the  first  geophone  to  the 
second.  The  speed  was  calculated  at  170  m/s. 

The  robot  was  kept  stationary,  and  the  geophones  were  rotated  around  it. 
Readings  were  taken  every  30°.  Measured  beam  pattern  was  calculated  by  the  following 
means:  the  maximum  received  voltage  at  each  bearing  was  noted  and  squared;  the 
squared  voltage  from  each  angle  was  then  plugged  into 

BeamPatterniO )  =  lOxlog lQ[signal(0) / signal(0° )] .  (5.1) 


This  beam  pattern  was  then  compared  to  a  theoretical  beam  pattern.  As  stated 
above,  the  robot  was  modeled  as  a  two-element  array  with  N=2  (number  of  elements). 
Spacing  between  elements  was  d=0.  27m.  The  wave  number  “k”  was  calculated  using: 

k=2tflc  (5.2) 


Where  f  =  frequency,  and  c  =  speed  of  sound.  The  directional  factor,  H,  was  calculated 
using: 


mo)  = 


1  sin[(  N  /  2)kd  sin  0] 
N  sin[(l/2)fa/sin#] 


(KFCS,  2000) 


(5.3) 


Where  N  =  number  of  elements,  k  =  wavenumber,  d  =  spacing  between  elements,  and  0  = 
bearing  of  geophones  relative  to  the  front  of  the  robot.  Theoretical  beam  pattern  was  then 
calculated  using: 

BeamPatterniO)  =  20 x  log10 (H(0))  (KFCS,  2000).  (5.4) 

The  results  are  shown  in  Figure  21.  Measured  beam  pattern  is  in  black,  and  theoretical 
beam  pattern  is  in  gray.  The  numbers  around  the  outside  of  the  plot  indicate  geophone 
bearing  relative  to  the  robot  in  [deg],  and  the  magnitude  scale  is  in  [dB]. 
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Figure  21.  Measured  vs.  Theoretical  beam  pattern  at  the  softball  field. 


The  theoretical  and  measured  beam  patterns  match  up  well  enough  to  not  disprove 
that  the  two-element  assumption  was  correct.  At  the  softball  field,  a  100  Hz  single  cycle 
tone  burst  was  the  only  signal  used.  Tests  using  different  signals  would  be  performed 
later  at  Del  Monte  Beach. 

It  was  desired  to  know  if  the  source  had  the  same  properties  in  different  media,  so 
further  beam  pattern  testing  was  performed  at  Del  Monte  Beach.  The  testing  was  done 
according  to  the  same  protocol  as  at  the  softball  field.  Ground  truthing  was  also 
performed.  An  example  is  shown  below.  A  lead  brick  was  dropped,  and  each  geophone 
picked  up  the  signal.  As  explained  above,  the  time  between  received  signals  was  divided 
by  distance  between  the  geophones.  An  example  plot  of  data  from  ground  truthing  is 
shown  in  Figure  22.  The  response  of  the  near  geophone  is  shown  in  black,  and  the 
response  of  the  far  geophone  is  shown  in  gray. 
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Figure  22.  Ground  truth  example. 


The  time  between  peaks  was  0.0166  seconds.  Given  a  2m  distance  between 
geophones,  the  speed  of  sound  in  the  ground  was  120m/s.  Further  trials  confirmed  that 
speed. 

Beam  patterns  were  also  calculated  the  same  way  at  the  beach  as  at  the  softball 
field.  Figure  23  shows  the  beam  pattern  for  a  single  -cycle  tone  burst  at  100  Hz.  The 
geophones  were  partially  buried. 
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Figure  23.  Beam  pattern  at  the  beach,  100  Hz,  6m  distance. 


Theory  and  measured  matched  up  almost  exactly.  With  that,  it  was  seen  that  the 
robot  acted  as  an  omni-directional  source.  Further  testing  involved  increasing  the  drive 
frequency.  This  was  done  to  see  if  the  robot  could  be  estimated  as  a  two-element  array. 
Figure  24  shows  the  beam  pattern  for  a  single-cycle  tone  burst  at  200  Hz.  Measured  is  in 
black,  and  theoretical  is  in  gray. 
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Figure  24.  Beach  beam  pattern  at  200  Hz,  6m  distance. 


Geophone  bearing  in  [deg]  relative  to  the  robot  is  shown  by  the  numbers  outside 
the  graph,  and  beam  pattern  in  [dB]  is  shown  on  the  magnitude  scale.  As  can  be  seen,  the 
measured  beam  pattern  does  not  match  up  to  the  theoretical  beam  pattern.  The  disparity 
between  measured  and  theoretical  grew  even  larger  as  frequency  was  increased  to  200  Hz 
and  250  Hz.  There  were  two  possible  conclusions  from  this:  a)  the  tone  burst  had  mixed 
P-waves  and  Rayleigh  waves  at  this  distance,  and  tests  should  be  repeated  at  a  greater 
range  which  would  allow  the  wave  types  to  separate;  or  b)  the  robot  could  not  be 
approximated  as  a  two-element  array  and  displayed  no  directionality  as  a  source.  This 
second  conclusion  would  be  most  unfortunate,  as  directionality  would  allow  beam¬ 
steering.  Beam  steering  would  greatly  increase  the  accuracy  of  mine  hunting  because  the 
operator  could  get  a  mine  bearing  as  well  as  distance. 
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2.  Wave  Classification 


The  next  task  was  to  determine  which  type  of  waves  propagated  from  the  robot. 
This  was  done  using  the  Stanford  Research  Systems  SR785  Signal  Analyzer. 
Measurements  were  taken  using  a  continuous  wave  frequency  scan.  This  meant  that  the 
SR785  put  out  a  continuous  wave  that  started  at  a  low  frequency  (ex.  50  Hz)  and  ended  at 
a  higher  frequency  (ex.  250  Hz).  The  frequency  scan  took  approximately  30  seconds  to 
complete.  Voltage  response  of  the  geophones  was  recorded  at  each  frequency.  The 
SR785  was  capable  of  recording  signal  magnitude  and  phase. 

After  doing  ground  truthing,  the  sound  speed  of  the  day  was  found  to  be  100  m/s. 
It  was  necessary  to  make  approximations.  Referencing  equation  (2.5),  Poisson’s  ratio  for 
sand  is  typically  v=  0.25  (Fitzpatrick,  1998).  This  gives  a  ratio  of  shear  wave  speed  to 
longitudinal  speed,  y  =  0.429.  Then,  solving  by  iteration,  it  was  found  that  the  ratio  of 
Rayleigh  wave  speed  to  the  shear  wave  speed,  k=0.898.  This  then  indicated  a  Rayleigh 
wave  speed  of  0.898* lOOm/s  =  90  m/s  =  Cr.  For  a  Poisson’s  Ratio  of  v=0.  25,  the  P- 
Wave  speed  is  approximately  3  times  the  Rayleigh  wave  speed  (Gaghan,  1998).  That  led 
to  a  P-Wave  speed  of  270  m/s. 

The  next  experiment  was  to  determine  if  Rayleigh  waves  were  excited.  To  do 
this,  radial  and  vertical  phases  were  recorded  as  a  function  of  frequency,  while  doing  a 
frequency  scan  with  the  SR785.  Both  the  radial  and  vertical  components  were  measured 
by  the  same  geophone.  The  set  up  for  the  phase  difference  experiment  is  shown  in  Figure 
25. 
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Figure  25.  Set  up  for  phase  difference  measurements. 
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Figures  26  a)  and  b)  show  phase  difference,  radial  component  minus  vertical  component, 
of  the  signal  put  out  by  the  robot.  The  geophone  was  oriented  so  that  radial  geophone 
motion  towards  the  robot  was  recorded  as  positive  radial  phase  measurement. 
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Figure  26  a)  and  b).  Radial  component  of  phase  minus  vertical  component  at  two  distances. 
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Frequency  [Hz]  is  on  the  x-axis,  and  phase  difference  [deg]  is  on  the  y-axis.  Figure  26  a) 
was  for  a  geophone  distance  of  4m  from  the  robot.  Figure  26  b)  was  for  a  geophone  12m 
from  the  source.  At  low  frequencies,  both  figures  show  a  center  on  -90°  phase  difference 
between  radial  and  vertical  signal.  A  -90°  difference  indicates  retrograde  elliptical 
motion.  Retrograde  elliptical  motion  indicates  the  predominance  of  Rayleigh  waves.  It 
is  interesting  that  the  maximum  frequency  at  which  Rayleigh  waves  predominate 
decreases  with  increasing  range.  At  the  higher  frequencies,  both  a)  and  b)  show 
significant  change,  with  both  actually  shifting  their  center  to  0°  .  This  shift  is  consistent 
with  the  phase  of  P-waves,  and  clearly  indicates  interference  from  P-waves.  Regardless 
of  the  interference  however,  Figure  26  shows  that  Rayleigh  waves  were  excited  by  the 
robot  source.  That  was  a  very  good  result  in  terms  of  mine  detection,  as  Rayleigh  waves 
were  the  preferred  means  of  detection. 

The  next  test  was  to  further  see  if  only  Rayleigh  waves  were  being  excited.  This 
involved  measuring  the  received  vertical  signal  at  various  distances  over  a  frequency 
scan.  Figure  27  shows  the  data  for  two  distances  14m  (black)  and  22m  (gray). 
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Figure  27.  Vertical  signal  magnitude  at  14m  and  22m  for  a  frequency  scan. 
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The  data  that  was  gathered  was  not  encouraging.  At  14m  there  were  nulls  every  4  Hz, 
and  at  22m  there  were  nulls  every  3  Hz.  This  led  to  the  suspicion  of  destructive 
interference  from  P-Waves  interfering  with  Rayleigh  waves.  Figure  28  represents  what 
was  suspected. 


Source  Receiver 


Figure  28.  Bottom  bounce  of  P-waves  leading  to  destructive  interference  at  the  receiver. 


Further  confirmation  of  the  destructive  interference  hypothesis  was  wanted.  For  a  P- 
wave  traveling  along  a  single  bottom  reflected  path,  the  phase  of  the  P-wave  is  given  by: 

Phase  =  -2k  p  Jh2+(-^)2  +  n  (5.5) 

Where  kp  is  the  wave  number  of  the  P-wave,  h  is  the  height  of  the  layer  the  wave  is 
traveling  in,  and  1  is  the  horizontal  distance  the  wave  travels. 

The  phase  difference  between  the  Rayleigh  and  P-waves  is  then: 

PhaseDiff  =  kRl  -  2 kp  Jh2+(^)2  +  n  (5.6) 

Where  kR  is  the  wave  number  of  the  Rayleigh  wave.  Destructive  interference  will  occur 
at  a  phase  difference  of  A=n,3n,5n,...  =  (2m+l);  m  =  0,  ±1,±2,....  After  substituting  in 

k  =  27if/c  and  some  algebraic  manipulation: 
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Where  Af  =  frequency  gap  between  nulls,  cR  =  Rayleigh  wave  speed,  and  cp  =  P-wave 
speed.  Layer  thickness  was  set  at  h  =  3m,  according  to  the  advice  of  Tim  Stanton,  an 
associate  researcher  in  the  Oceanography  Department  at  NPS.  Calculated  values  for 
frequency  gaps  between  nulls  were  10  Hz  at  14m  and  6  Hz  at  22m.  These  values  did  not 
exactly  agree  with  estimates  from  looking  at  the  data.  However,  they  were  in  the  same 
ballpark  as  the  measured  values.  Also,  the  calculated  theoretical  results  did  show  the 
same  decrease  in  frequency  between  nulls  as  distance  increased  that  the  data  showed. 
The  results  of  the  calculations  were  strong  enough  to  conclude  that  destructive 
interference  between  P-waves  and  Rayleigh  waves  was  occurring.  Destructive 
interference  could  decrease  Rayleigh  wave  amplitude,  and  could  also  lead  to  false  alarms 
due  to  P-wave  detection  by  sensors.  Destructive  interference  and  P-wave 
propagation/bottom  bounce  therefore  would  seriously  hinder  any  attempts  at  target 
detection. 


C.  TARGET  DETECTION 


1.  Detection  with  Oscilloscope 


Despite  the  discouraging  results  from  the  wave  classification  measurements,  a 
target  search  was  tried  just  in  case  it  actually  worked.  Figure  29  shows  the  set-up  used. 
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Figure  29.  Set-up  for  target  search  experiment. 


For  the  first  trial,  the  geophone  was  2m  from  the  robot,  and  the  target  was  4m  from  the 
geophone  (6m  from  the  robot).  However,  many  different  distances  between  the  robot  and 
the  geophone,  and  the  geophone  and  the  target  were  sampled.  The  target  was  a  large  air 
cylinder,  which  had  a  mass  of  55  kg.  It  was  buried  just  beneath  the  surface  of  the  sand. 

It  was  hoped  that  the  geophone  would  first  detect  the  initial  signal  from  the  robot,  and 
then  detect  an  echo  signal  from  the  target. 

Data  from  the  target  search  experiment  was  processed  using  “vector  polarization 
filtering”.  Vector  polarization  filtering  is  a  signal-processing  technique  for  suppressing 
all  but  the  Rayleigh  waves  in  the  received  signal  (Sheetz,  2000).  It  is  based  upon  the  fact 
that  for  a  Rayleigh  wave,  the  ground  motion  is  elliptical  and  so  the  radial  and  vertical 
components  are  90°  out  of  phase.  For  a  P-wave,  the  radial  and  vertical  components  are  in 

phase.  If  Vr  (!)  and  Vv  (t)  are  the  phasor  representations  of  the  radial  and  vertical  ground 
motion  signals,  respectively,  then  the  following  complex  quantity,  called  the  “complex 
power”,  is  formed: 

* 

Pr,y(t)  =  Vr(t)»Vy(t)  .  (5.8) 

Where  Pr  v  (t)  is  the  complex  power,  and  V  is  the  complex  conjugate  of  the 

radial  ground  motion  signal.  The  imaginary  component  of  the  complex  power  captures 
the  entire  Rayleigh  wave  signal  and  tends  to  suppress  the  P-wave  signal.  The  (complex) 
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phasor  representations  of  the  radial  and  vertical  ground  motion  signals  are  formed  from 
the  real-valued  collected  data  by  a  Hilbert  transform  operation  using  a  MATLAB  M-File. 
The  code  for  the  program  is  shown  in  Appendix  D 

Figure  30  below  is  for  a  set-up  with  no  target.  Figure  31  is  for  a  set-up  with  a 
target  in  place,  6m  from  the  robot.  LT  Sheetz  had  success  with  target  search  at  the  4-6m 
range  for  target  search,  so  that  is  why  6m  was  chosen.  For  both  tests,  the  drive  was  a 
single-cycle  tone  burst  at  100  Hz,  and  drive  voltage  was  47  V  peak-to-peak.  Readings 
were  recorded  on  the  Hewlett  Packard  Infinium  oscilloscope,  with  64  averaging. 
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Figure  30.  Imaginary  part  of  received  power  at  2m  from  source,  no  target. 
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Figure  31.  Imaginary  part  of  received  power  at  2m  from  source,  with  a  target  at  6m. 


For  a  Rayleigh  wave  speed  of  90m/s,  the  initial  signal  was  supposed  to  hit  the 
geophone  at  time  t=0.044  seconds.  Both  plots  show  the  initial  signal  arriving  at  that 
expected  time.  On  the  “with  target”  experiment,  a  return  signal  was  expected  at  time  t  = 
0.110s.  Fiowever,  the  plot  for  the  “no-target”  experiment  is  nearly  identical  to  the  plot 
for  the  “with-targef  ’  experiment,  even  out  at  time  =  0.1 10  seconds.  This  means  that  there 
was  no  target  detection.  Experiments  at  different  distances  showed  the  same  negative 
results. 


2.  Detection  with  Human  Ear 


For  curiosity  sake,  detection  by  listening  was  attempted.  The  geophone  panel  was 
used  to  transmit  a  geophone’s  received  signal  into  a  stereo  headphone  amplifier. 
Headphones  were  then  plugged  into  the  amplifier,  and  the  operator  could  actually  hear 
the  received  “thump”  at  the  geophone.  Thumping  could  still  be  heard  clearly  at  ranges  of 
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20+  meters  away  from  the  source.  However,  there  was  a  rising  tone  as  the  geophone  got 
further  away  from  the  source.  This  indicated  that  the  lower  frequencies  were  reaching 
the  geophone  before  the  higher  frequencies,  i.e.,  the  waves  arriving  at  the  geophone  were 
dispersive.  This  dispersion  is  probably  due  to  and  increasing  wave  speed  with  the 
increasing  amount  of  compaction,  rigidity,  and  moisture  content  with  depth  in  the  sand. 

Target  detection  by  this  method  was  then  attempted.  It  was  hoped  that  two 
separate  signals  could  be  heard  through  the  headphones:  the  initial  signal  from  the  robot; 
and  the  echo  signal  from  the  target.  Attempts  were  unsuccessful  as  only  one  thump  was 
heard. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Landmines  continue  to  be  a  threat  to  both  military  and  civilian  communities 
throughout  the  world.  Current  methods  of  detection,  while  better  than  nothing,  could 
certainly  be  improved.  Seismic  SONAR  is  a  promising  new  technology  that  may  help 
save  countless  lives. 

The  goal  of  this  thesis  was  to  advance  Seismic  SONAR  development  by 
introducing  a  mobile  source  that  could  be  easily  used  in  practical  applications.  Used  for  a 
mobile  source  was  a  small  tracked  vehicle  with  dual  inertial  mass  shakers  mounted  on 
top.  The  source  accurately  transmitted  the  shaker  signal  into  the  ground,  and  its  mobility 
made  it  a  practical  choice  for  field  operations.  It  excited  Rayleigh  waves,  as  desired,  but 
also  generated  undesirable  P-waves  and  was  not  directional  at  all.  It  proved  incapable  of 
finding  a  target. 

There  were  three  major  reasons  why  this  particular  source  failed  for  mine 
detection.  First,  it  must  be  possible  to  selectively  generate  Rayleigh  waves,  while 
suppressing  P-waves.  Interference  from  bottom-bouncing  P-waves  be  destructive  and 
greatly  confused  signal  processing.  Second,  more  aperture  is  needed  so  that  a  steer-able, 
directional  beam  can  be  produced.  This  would  allow  the  operator  to  know  range  and 
bearing  for  buried  mines.  The  directionality  problem  may  be  solved  by  simply 
preferentially  exciting  Rayleigh  waves  and  excluding  P-waves.  Third,  a  stronger  source 
is  needed.  This  could  be  accomplished  by  either  a)  using  more  shakers,  or  b)  adding 
more  weight  to  the  robot.  This  source  could  be  made  to  be  effective  if  mounted  with  an 
array  of  shakers.  CPT  Sheetz  was  able  to  produce  a  highly  directional  beam  with  a 
strong  enough  signal  for  target  detection  and  bearing  by  using  a  7  element  linear  array  of 
dual  Aura  Bass  Shakers.  This  same  concept  could  be  applied  to  tracked  vehicles.  Figure 
32  shows  a  possible  design  for  mounting  an  array  of  shakers  on  a  tracked  vehicle. 
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Figure  32.  Proposed  design  for  further  testing. 


The  array  could  easily  be  lifted  and  lowered  by  means  of  mechanical  arm,  and 
mobility  of  the  robot  would  not  be  compromised.  Each  shaker  could  be  at  a  45°  angle 
with  respect  to  the  ground.  This  would  preferentially  excite  elliptical  particle  motion, 
thus  selectively  producing  Rayleigh  waves  and  excluding  P-waves.  The  concept  of  a 
tilted  source  selectively  producing  Rayleigh  waves  has  already  been  shown  by  LT 
Gaghan  and  can  work. 

Seismic  SONAR  is  still  a  promising  means  of  mine  detection,  and  it  has  been 
shown  to  work  with  a  fixed,  buried  array  of  sources.  Applying  the  concept  of  a  large 
linear  array  and  combining  it  with  a  mobile  vehicle  should  bring  seismic  SONAR  into  the 
field  of  military  use.  This  thesis,  while  not  actually  achieving  a  usable  source,  showed 
that  a  tracked  vehicle,  with  improvements,  is  a  means  of  mine  detection  with  high 
potential. 
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APPENDIX  A.  MANUFACTURER  SPECIFICATIONS  FOR  AURA 

BASS  SHAKERS 


Mir.  Sugg.  Retail  $199.95 


WAS 

w 
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Bass  Shakers”  Let  You 
Feel  All  the  Bass  without 
Breaking  the  Sound  Barrier 

The  3ass  Shaker  is  a  t'ansducer  that  gereratss  the 
sensation  cl  sound  by  vibration,  no:  by  moving  air  The  result  is  a  big  bass  effect  wtinoul 
a  high  cressure  revel  hat  could  distort  sound  or  blow  your  speakers  They  are  also  great 
lor  adding  bass  tc  vehicles  with  poor  sourd  insulation,  sucn  as  trucks.  Mount  these  Bass 
Shavers  unoer  your  iron!  seats  ard  you  will  teel  the  punch  of  he  drums  anc  the  kick  o'  the 
percussion,  as  it  you  were  on  stage  or  in  tha  stucio  when  the  music  was  recorced1  Hooks 
up  easily  to  any  standard  audio  ampli'ier  output  channel,  and  can  oe  used  with  existing 
suowopiers  for  an  extra  kick.  Each  Bass  Shaker  requires  only  25W  RMS  of  oower 
Mode!  AST1B4  Item  No.  B-40070-492383  S/H 


Specifications: 

Frame  Size: . 

Magnet  Type: . 

Power  Rating: . 

Force.  Nominal: . 

Weight: . 

Resonance  Frequency  (fo): 
Usable  Frequency  Range: . 


AST-1B-4  Bass  Shaker 
...  4.75"  dia.X  2.35"  h 
....  Ceramic 
...  Nominal  25  W 
....  10lbf(44.5  N) 

_ 2.5  lbs  each 

42  Hz 

....20-100  Hz 
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APPENDIX  B.  MANUFACTURERS  SPECIFICATIONS  FOR 

GEOPHONES 


SM-11  Geophone 


-  geophone  high  BJSijrtaua,  aver  Sflfi  Hz 

providing  wirt«  bpmjwmtn  dKa  nullable-  Tar  lp  fa 
1  mi  data  sampling 

■  Can  be  apnratn::  In  any  O'lnntallijfi 

■  liflh  culpul  Hvrogn  rhu  UE4  0*  a  ipSCml  Magnet 
and  czar  design 

■  Ruauzd  mechanical  Darainjclio*i  can  wnrstantl 
sen**  shocks 

■  J-ySB#  Heihcd  warranty 


The  Shfl-I  I  gbnphnna  it  nbitatee  I'flr  use  m  arlchded  'rHEUCTicy  higlT-re5.;;iLmc«l  Surveys.  It  has  a 
na1cra:  IreqLei-ry  nl  30  Hr  snd  a  ecuIXH/S  rvtrjcncy  nf  over  SM  Hr  fiirivir  -  *  kehsea  sL'Iatie  lor 
new  v«llh  l-his  lamping  izcnrdmg  syslams  Iha  ar  i  kpeda  mapnellE  oiruc  nan#t  1h#  ouftiU! 
al  Ihii  gcnphqnfl  nigh*  H*ih  harmol  HO-Hr  gecptwrrB*  imsur  hS  Wequnle  signal  st'anglh.  The  li|gn 
nadura  IrnquonCf  spnng  design  jh;c-  mims  Ibis  geoptinnrs  rg  W  OHd  n  flay  Unenlalisn  tyr-tical. 
hanzzmtal  or  ivftTJiiJ, 

The  £M-vi  am  bn  mslBUfld  in  the  WO  Sehttjr  PE-df  Lind  case 
T,pital  a  L|L'i-:jlic-n  h>gh TfeMlijOan  wjhBmijQraivi  i*elflec*or  aluaiei 
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APPENDIX  C.  MANUFACTURERS  SPECIFICATIONS  FOR 

ACCELEROMETERS 


'  Calibration  Data  Card 
SHEAR  ACCELEROMETER 

MODEL  *  _ 

SERIAL  *  59 C  4 


^titCallbralion  Data  Card 

rrtvyf  SHEAR  ACCELEROMETER 
MODEL  a  73S'?  BO  i 

SERIAL  *  ? 


/OLTAGE  SENSITIVITY  fO.OZ 
-REQOENCY  RANGE:  /  7Y?  H, 

x/TPLrr  bias  i  fvfi  ■  S'-  V _ v 

lain  ‘2-2f  ^V  fW 

For  further  information  please  refer  to 
Calibration  Certiteale 

PfB 

Nf«i  »  I 

3425  Walden  Ave  •  Depew,  NY-14043 
(716)664-0001 


VOLTAGE  SENSITIVITY:  tOi-V _ mVilg 

TREOUENCY  RANGE  jlO  7000  u- 

OUTFLIT  BIAS  LEVEL:,  g  / _ v 

Date:  f*7l"VY  gy  C  t/ , 

For  farther  mtorntalion  please  refer  to 
CaUxaton  oeitmcate. 

PIB 

■v-wm 

3425  Walden  Ave  •  Oepew.  NY- 14043 
(716)684-0001 
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APPENDIX  D.  MATLAB  CODE  FOR  “HILBRT.M” 


%  Name: 

Hilbrt . m 

%  Author: 

LT  Mike  Fitzpatrick 

%  Updated:  05/30/03  by  ENS  Doug  MacLean 

%  Description: 
Hilbert  analysis  of  the 

This  subroutine 

conducts 

a 

o 

o 

polarization  filtering. 

geophone  signals 

for 

vector 

o 

o 

output  the 

This  is  a  versatile  program 

meant 

to 

o 

o 

imaginary  and  real 

parts 

of 

a 

received  signal .  There 


are  various  user-selectable  features 


%%%Input  Parameters%%% 

Range=l;  %Turns-on  plotting  with 

range  axis 


wavespd=296; 

%Wavespeed  [ 

ft/s] 

%%%Date,  Directory, 

&  File%%% 

clc,disp( ' ***Hilbert 

Analysis 

Subroutine***  '  ) 

;  di sp  (  '  '  )  ; 

if  exist (' date ') ~=1 , 

dir  *  .  , 

date=input  (  ' Enter 

the  date 

— >  ' ,  '  s  1 )  ; 

end 

clc,disp( '***Hilbert 

Analysis 

Subroutine***  '  ) 

,  dir  *. 

directory=input  (' Enter  the  di 

rectory  — >  '  ,  ' 

s'  )  ; 

cd  (directory) ; 

clc,disp( ' ***Hilbert 

Analysis 

Subroutine***  '  ) 

,  dir  *.mat; 

file=input (' Enter  filename  — 

>  '  ,  ' s ' ) ;  disp ( 

'  '); 

data=xlsread (file)  ; 

transform=hilbert (data)  ; 
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Pwr ( : , 1 ) =con j (transform (:,!))  . * transform ( : ,  2 ) ; 


directory=input (' Enter  the  directory  for  time  vector  --> 

'  ,  '  s  '  )  ; 

cd  (directory) ; 

clc, disp ( ' ***Hilbert  Analysis  Subroutine***'),  dir  *.mat; 

timeFile=input (' Enter  filename  for  time  vector— >  ','s'); 

disp  (  '  '  )  ; 

t=xlsread (timeFile) ; 

%%%Plotting%%% 
figure  (1) 
plot (t, real (Pwr) ) 
title ('Real  Power  -->  ') 
ylabel ( 'Amplitude ' )  ,  grid 

oooooooooooooo 

figure  (2) 
plot  (t, imag (Pwr)  ) 
title  (  ' Imaginary  Power ' ) 
xlabel('Time  [ s ] ' ) 
ylabel ( 'Amplitude ' )  ,  grid 
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